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Tethered  Satellite  System  Stability 

D.-C.  Liaw  and  E.H.  Abed 

Department  of  Electrical  Engineering 
and  the  Systems  Research  Center 
University  of  Maryland 
College  Park,  MD  20742  USA 

Abstract 

Issues  of  stability  of  the  Tethered  Satellite  System  (TSS)  during  station-keeping,  de¬ 
ployment  and  retrieval  are  considered.  The  basic  nonlinear  equations  of  motion  of  the 
TSS  are  derived  using  the  system  Lagrangian.  Using  the  Hopf  bifurcation  theorem,  ten¬ 
sion  control  laws  are  established  which  guarantee  the  stability  of  the  system  during  the 
station-keeping  mode.  A  constant  angle  control  method  is  hypothesized  for  subsatellite 
deployment  and  retrieval.  It  is  proved  that  this  control  law  results  in  stable  deployment 
but  unstable  retrieval.  An  enhanced  control  law  for  deployment  is  also  proposed,  which 
entails  use  of  the  constant  angle  method  followed  by  a  station-keeping  control  law  once 
the  tether  length  is  sufficiently  near  the  desired  value.  Simulations  are  given  to  illustrate 
the  conclusions. 

I.  Introduction 

The  Tethered  Satellite  System  (TSS)  [l]-[8]  consists  of  a  satellite  and  subsatellite 
connected  by  a  tether,  in  orbit  around  the  Earth.  Many  potential  applications  of  the 
TSS  have  been  proposed,  including  deployment  of  scientific  instruments  and  study  of  the 
Earth’s  magnetic  field  [1],  [12].  Station-keeping,  deployment  and  retrieval  of  payloads  are 
the  three  major  modes  of  operation. 

Arnold  [2]  proposed  a  constant  angle  method  for  deployment  and  retrieval  of  subsatel¬ 
lite  of  the  tethered  satellite  system.  In  [2],  the  satellite  and  subsatellite  are  modeled  as 
point  masses  and  the  tether  is  assumed  massless  and  of  length  small  compared  with  the  ra¬ 
dius  of  the  satellite’s  orbit.  Based  on  these  assumptions,  Arnold  obtained  an  approximate 
model  of  the  TSS  by  applying  the  gravity-gradient  method  and  argued  that  this  constant 
angle  scheme  would  result  in  stable  deployment  and  unstable  retrieval. 


One  goal  of  this  paper  is  to  give  a  proof  of  the  validity  of  these  conclusions.  First,  how¬ 
ever,  general  dynamic  equations  for  the  TSS  are  derived  by  using  the  system  Lagrangian. 
The  applied  tension  control  force  is  assumed  to  be  the  only  external  force  acting  on  the 
TSS.  Next  we  consider  stabilization  of  the  TSS  during  the  station-keeping  mode,  in  which 
the  tether’s  length  is  regulated  to  remain  nearly  fixed.  We  observe  from  linear  analysis 
at  the  system  equilibria  the  presence  of  two  pairs  of  pure  imaginary  eigenvalues  in  the 
absence  of  feedback.  This  suggests  the  possibility  of  librations  superimposed  upon  the 
orbital  motion.  It  is  found  that  nonlinear  stability  analysis  is  needed  to  study  stability 
and  stabilization  of  the  TSS  during  station-keeping.  The  program  of  [15],  which  considers 
Hopf  bifurcation  control  algorithms,  is  employed  to  derive  stabilizing  control  laws.  Both 
linear  and  nonlinear  feedback  controls  are  achieved  which  guarantee  asymptotic  stability. 

Viewing  the  tether  length  as  an  input  variable  for  deployment  and  retrieval  of  the 
TSS,  a  constant  in-plane  angle  control  scheme  is  considered  next.  Within  this  setting, 
we  prove  stability  of  constant-angle  deployment  and  instability  of  constant-angle  retrieval. 
This  is  achieved  through  the  construction  of  appropriate  Liapunov-like  functions  and  by 
appealing  to  the  finite-time  stability  theory.  A  new  control  strategy  for  deployment  of 
the  subsatellite  is  also  proposed.  This  control  law  consists  of  the  constant  angle  scheme 
followed  by  the  stabilizing  station-keeping  control. 

Finally,  simulation  results  are  given  to  demonstrate  the  analytical  conclusions  of  the 
paper. 

Notation 

E  -  Earth 
S  -  Satellite 

m  -  Subsatellite  and  its  mass 
G  -  Gravitational  constant 

M,  ms  -  Mass  of  the  Earth,  mass  of  the  satellite 

(ajm,ym,  zm)  -  Earth-based  rotating  Cartesian  coordinates  of  subsatellite,  with  zm  in  the 
local  outgoing  vertical  direction,  and  xm  in  the  direction  of  motion  of  the  satellite 
in  its  orbit  (see  Figure  1) 
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(£m,ym,£m)  -  Inertial  coordinates  of  subsatellite 
(xa,y3,zs)  -  Inertial  coordinates  of  the  satellite 

-  Constant  angular  velocity  of  the  satellite  in  circular  orbit 
0,  <f>  -  In-plane  angle  and  out-of-plane  angle  of  subsatellite  relative  to  local  vertical 
:=  <j>,  u?e  :—$,£-  Tether  length,  v  :~  £ 
r0,  rm  -  Radius  of  the  satellite  orbit,  radius  of  subsatellite  orbit 
(•)*  -  Evaluation  at  equilibrium 
T0,  -  Torques  in  directions  #,  <f> 

Ft  -  Force  along  tether  (i.e.,  the  applied  tension  force) 

F  :=  +  -y  +  Fit ,  where  a  hat  indicates  a  unit  vector  in  the  given  direction 

II.  System  Dynamic  Equations 

The  coordinate  system  of  a  typical  tethered  satellite  system  is  depicted  in  Figure  1. 
Referring  to  this  coordinate  system,  we  make  the  following  simplifying  assumptions:  The 
satellite  and  the  subsatellite  are  point  masses.  Their  masses  are  related  as  ms  m,  and 
hence  the  center  of  mass  of  the  TSS  may  be  taken  to  coincide  with  the  satellite.  The 
tether  is  massless  and  rigid,  and  the  gravitational  attraction  between  the  subsatellite  and 
the  satellite  is  neglected.  Finally,  the  TSS  experiences  no  aerodynamic  drag  forces,  and 
the  satellite  is  in  a  circular  orbit  around  the  Earth. 

y 


Figure  1.  Coordinate  System  of  the  TSS 
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It  is  obvious  from  Figure  1  that  we  have  following  relationships  (see  the  notation  list 
in  Section  I). 


xm  —  i  cos  <f>  sin  9 
Vm  =  i  sin  (j) 

zm  —  ^0  +  £  COS  <f>  COS  9 

rm=r o  +  +  2ro£  cos  (f>  cos  9 

The  inertial  coordinates  of  the  subsatellite  and  the  satellite  are  given  as 


xr 

Vm 

Zm 


cos  Of  0  sin  ilt 
0  1  0 
—  sin  f It  0  cos  f It 


Xm 

Dm 

Zm. 


(1) 

(2) 

(3) 

(4) 


(5a) 


1  xs  > 

/ 

Vs 

\  Zs  / 

\ 

cos  Sit  0  sin  Slt\  /  0 

o  i  olio 

—  sin  Sit  0  cos  f  It  J  \r  0 


(56) 


where  an  implicit  choice  of  time  reference  is  understood. 

The  kinetic  energy  of  the  TSS  is 

T._.  1  x2s  1  .  x2  ^2  . 

RE=  +  ya  +  Za)  +  -m(xm  +  ym  +  zm). 


(6) 


Using  Eq.  (5),  this  may  be  rewritten  as 

I<E  =  \msSl2r2  +  \m{£2  +  +  l2  cos2  <j>(9  +  SI)2  +  Sl2r 2 

Z  Z 

+  2firo^cos  (j)  sin#  —  20rot?sin<^>sin#<^  +  2Slro£  cos  cf>  cos  9(9  +  fl)}. 


(7) 


The  potential  energy  of  the  TSS  is  due  to  gravity  and  is  given  by 

GMrris  GMm 


PE  = 


r  o 


(8) 


The  satellite  is  in  a  circular  orbit,  implying  that  it  is  in  a  zero-g  orbit.  Hence 

GMrris  o2 
- 5 —  =  m3U  r0, 


(9) 
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i.e.,  we  can  take  GM  =  ftVg. 

By  writing  the  system  Lagrangian  L  —  KE  —  PE  and  invoking  the  Lagrangian  for¬ 
mulation  of  the  system’s  dynamics,  the  equations  of  motion  of  the  TSS  are  obtained  as 

j>  .  . 

Te  =  mi 2  cos2  4>{9  +  2-(0  +  ft)  —  2tan <f>(9  +  f V)<f> 

ft2r0  sinff  _  rjL. 

I  cos  <f> 

l  . 

=  m£2{^>  +  2-<^  +  cos  <j>  sin  <j)(9  +  ft)2 


02r-A  r3 

4 - —  cos0sin<£(l - |-)}  (11) 

i  rm 

=  m{i!  —  £(<f>)2  —  £cos2  <f>(9  +  ft)2 

OW  3 

H - ^ - ft2r0  cos  ^  cos  #(1 - |-)}  (12) 


For  the  case  r o  £,  we  have  rm  ~  r0.  Moreover,  in  this  case  Eq.  (4)  implies 

r3  £ 

1 - ~  3  cos  <^>  cos  9 — .  (13) 

r3m  ^0 

Hence,  the  approximate  equation  of  motion  for  the  system  for  the  case  r0  £  is  obtained 
as 

F  =  m£{£  -  -  £  cos2  <f>(9  +  ft)2  +  m2  -  3ft2 1?  cos2  <f>  cos2  9} 

+  m9{9£  cos  <f>  +  2(9  +  ft)(f  cos  <f>  —  £(j>  sin  <j>)  +  3 £G?  cos  9  cos  <f)  sin  9} 

+  m(j>{£<i>  +  2 £(j)  +  £  cos  <f>  sin  <f>(9  +  ft)2  +  3£ft2  cos2  9  cos  <j)  sin  <j)}. 


This  agrees  with  a  result  of  Arnold  [2].  Note,  however,  that  we  do  not  require  r0  >  £ 
below. 

III.  Stabilization  of  System  During  Station-Keeping  Mode 

To  facilitate  design  of  stabilizing  control  laws  for  the  TSS,  we  first  rewrite  the  system  of 
second  order  equations  (10)-(12)  in  state  space  form.  This  is  possible  under  the  assumption 
cos  (j>  ^  0  (i.e.,  <f)  /  ±f ).  In  addition,  we  assume  that  the  only  external  force  acting  on  the 
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TSS  is  the  applied  tension  Ft.  Setting  re  =  t#  =  0  in  (10)-(12)  and  denoting  T  := 
equations  of  motion  in  state  space  form  are  found  to  be 

(j)  —  LVtj,  -  - 

sin(2<f>l(oM -MZT - 

£  v  2 


Ci+  =  -  y  i  sin(2<^)(a;e  +  ft)2  -  ^  cos  0  sin  <X  1  “  y) 


0  =u>8 

2v  \  ft2rosin#.  rjj  . 

d><?  =  — ~f(u>6  +  ft)  +  2  tan  <f>(u)$  +  £l)u)4>  £  Cos  <f>  '  _ 


£  =  v 

Cl2r3£ 

v  —  £w\  +  £  cos2  <p{ui8  +  ft)2 - j*— 


Ft,  the 

(14) 

(15) 

(16) 

(17) 

(18) 


7* 

+  fl2ro  cos  0  cos  (f>{\ - 1-)  + 


rrn  171 


(19) 


With  the  tether  length  held  constant  {£  =  £*),  this  system  has  precisely  two  equilibrium 
points  (0,0,0,0,f  ,0)  and  (0,0, ^,0/*  ,0).  Moreover,  the  linearized  system  of  Eqs.  (14)-(17) 
at  such  an  equilibrium  is  found  to  have  two  pairs  of  pure  imaginary  eigenvalues,  given  by 


Ai, 


2  =  iiftyl 


+  S(1-3M. 


«-■ 


m,  0 


(20) 


-^3,4 


(21) 


where  rm,0  =  r0  +£*  at  9  =  0  and  rm,o  =  r0 -P  at  9  =  ir.  The  pairs  Ai>2  and  A3,4  of  eigen¬ 
values  are  associated  with  the  out-of-plane  and  in-plane  dynamics,  respectively,  suggesting 
that  the  TSS  may  undergo  librations  with  two  different  frequencies  superimposed  upon 
the  orbital  motion,  near  either  of  the  equilibria.  In  addition  [9],  the  set  <j>  =  0,w*  =  0  is  an 
invariant  manifold  for  Eqs.  (14)-(19),  regardless  of  the  tension  control  force  T.  Although 
this  implies  the  uncontrollability  of  the  system,  a  tension  control  law  still  can  be  designed 


to  guarantee  asymptotic  stability. 

Stabilizability  will  first  be  studied  at  the  equilibrium  point  x0  =  (0,0,0,0,£*,0)T.  Apply 
a  linear  tension  control  law  T  =  -m(U  +  M  +  +  M  +  hv),  where  kiy  *  =  1,  •  •  • ,  6 
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(22) 


are  constant  control  gains,  £:—£  —  £*  and 

(3rgl»  +  3r0l*2  +  l*3)ft2 
(r0  +  £*)2 

The  characteristic  equation  of  the  linearized  closed-loop  model  of  system  (14)-(19)  is  found, 
after  some  manipulation,  to  be 

(A2  +  a2)  (A4  +  &4 A3  -f-  61A2  +  62A  +  63)  =  0,  (23) 


where 


,4rg+6rj£*+4r0r2  +  r\0.5„ 

“■  = ( - WTff - } 

,  ,  (2.„r2  +  «*3)a2  2 sik2  , 

bi=h-  (ro+<.) 3 - F“-  +  4!i  ■ 

t4(3r2  +  3rJ**  +  r„<*2  )!22  2!!fc, 

2_  (rc+<*)>  “  <*  ’ 

(3rjj  +  Zrlt  +  r0r2)ii2  ,  (3rg  +  3rgl*  +  3rpr2  +  C3)n2 
3  (rii  f  i*  )2  1  3  (r0+£*)3 


(24) 

(25) 

(26) 
(27) 


By  applying  the  Routh-Hurwitz  test  to  Eq.  (23),  we  obtain  the  following  preliminary 
result. 

Lemma  1.  If  the  tension  control  force  is  given  as  T  =  m(—U  —  k\6  —  k20Je  —  k^£  —  k±v), 
then  the  linearized  closed-loop  system  of  Eqs.  (14)-(19)  at  the  equilibrium  point  x0  has  a 
pair  of  pure  imaginary  eigenvalues  and  four  eigenvalues  with  negative  real  parts,  provided 
that  the  gains  ki,i  =  1, ...  ,4  satisfy  the  following  two  conditions: 

(i)  h,b2,b3,k4  >  0 

(ii)  &4&1&2  —b\  —  k\bz  >0 

Here,  the  bi,i  —  1,2,3  are  as  given  in  (25)-(27). 

Note  that,  by  Eq.  (23),  the  system  (14)-(19)  has  an  uncontrollable  pair  of  pure 
imaginary  eigenvalues  ±i«i,  unaffected  by  linear  state  feedback..  It  is  easily  checked  that 
this  holds  even  with  feedback  of  states  (j>  and  The  stability  of  the  closed-loop  system 
cannot,  therefore,  be  determined  from  the  linearized  model.  That  is,  this  is  an  example  of 
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a  critical  case  in  nonlinear  stability.  Stability  results  for  Hopf  bifurcation  of  one-parameter 
families  of  nonlinear  systems  can  be  employed  in  studying  the  stability  of  critical  systems 
with  a  single  pair  of  pure  imaginary  eigenvalues.  For  a  discussion,  see  [15].  Two  types 
of  stabilizing  tension  control  laws  are  given  in  the  next  two  theorems  for  the  station¬ 
keeping  application.  Theorem  1  involves  the  use  of  linear  state  feedbacks  and  Theorem  2 
accounts  for  a  class  of  nonlinear  feedbacks.  Details,  proofs  and  simpler  results  in  the  form 
of  corollaries  can  be  found  in  [13]. 

Theorem  1.  If  a  linear  state  feedback  controller  T  of  the  type  specified  in  Lemma  1  is 
applied,  with 

(~2~  +  2  1  — 4" (&3  —  «3))di  —  -k4ail*d,2  <  0, 

then  the  equilibrium  point  rco  is  rendered  asymptotically  stable  for  the  system  (14)-(19), 
where  a\  is  as  given  in  (24),  and  where 


d\  =  ki  +  ~2^(,a2  +  4ai), 

d2  =  2a\(k2  ~  2£*£l)  H - ^  ^  y — — («3  +  4a^  —  k%), 

_  (3r3  +3rgr  +  r0l*2)ft2 

(r0+^*)3 

_  (3rp  +  3 r2l*  +  3 r0l*2  +  l*3)I22 
a3  (r0+^*)3 

(8r^*  +  14r3£*2  +  16r2£*3  +  9  r0t4  +  2£*5)tt2 
ai~  2(r0  + 1*)* 


(28) 

(29) 

(30) 

(31) 

(32) 


Theorem  2.  If  the  applied  tension  control  force  is  given  as  T  —  m(—U  —  k\6  —  k2u>o  — 
k^l  —  k±v  —  qi<t>2  —  I2 4>u<t>  —  430^)  with  ki ,  i  —  1,  . . .,  4  satisfying  the  conditions  of  Lemma 
1,  then  the  system  (14)-(19)  is  rendered  asymptotically  stable  if  the  condition 

/  4l  T  ( I"  “f"  43  \2  / 1  \  \  t  1  /  0*  U  \  A  ^  n 

( +  ( — - — K  -  (h  -  a3)—)d  1  +  y  (-42  -  l  h)d2  <  0 

holds.  Here,  di,  d2  and  ai,i  =  1, . . .  ,4  are  as  defined  above. 
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Similar  results  can  also  be  obtained  for  stabilizing  the  system  at  the  equilibrium  point 
(0, 0, 7r, 0, €*, 0)  [16].  The  details  are  omitted. 

IV.  Constant  In-Plane  Angle  Control 

In  this  section,  we  consider  deployment  and  retrieval  of  the  subsatellite  in  the  tethered 
satellite  system.  Viewing  i  as  an  external  control  input,  the  state  equations  (10)-(12)  of 


the  system  are 

0  =  ue  (33) 

2 i.  fiVosin#  .  rij  .  „ 

=  — — (ce#  +  f7)  +  2  tan  (f>( loq  +  - — - — — (1 - — )  (34) 

v  K,  COS  (p  ^ 7YI 

j>  =  w,j>  (35) 

2£  1  fi2rn  r3 

u<t>  =  g  sin(2<^)(^  +  ^)2 - j~  cos  0  sin  0(1  -  ^f-).  (36) 


At  an  equilibrium  point  of  (33)-(36),  if  one  exists,  we  have  Uq  =  —  0,  and 

t  is  given  (from  Eq.  (34))  by 


i  ro  \ 

2 cos 0*  Sin  (r^(0)3 


and  0*  must  satisfy  either 


sin  0*  —  0,  or 
ro 


COS  0* 


(1 


(rm(^)) 


Vwr)cos0*> 


(37) 


(38a) 

(386) 


where 

rrn(^)  :=  (r 0  +  ^2  +  2r0£cosO*  COS0*)1/2.  (39) 


Remark  1.  In  fact,  only  the  case  sin  0*  =  0  is  realistic.  To  see  this,  consider  momentarily 
the  possibility  (386),  which,  using  (33),  would  imply  that  at  equilibrium  i  obeys 


e  = 


(40) 
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Since  —  j  <  (/>*  <  f,  we  have  cos <f>*  >  0  (see  Figure  1).  Considering  the  possibilities 
0  <  <£*  <  f  and  —  f  <<(>*<  0  separately,  and  referring  to  Figure  1  for  the  relative 
magnitudes  of and  r0,  we  find  that  the  left  and  right  sides  of  (386)  are  then  of 
opposite  sign  unless  they  both  vanish.  Hence,  we  obtain  4>*  =  9*  =  ±-|,  implying  £  of  (40) 
would  be  infinite. 

□ 

In  the  light  of  Remark  1,  we  let  4 >*  =  0.  Eq.  (37)  now  implies  that  at  equilibrium  l 
satisfies 


i  = 


Qr0 


(i 


-)  sin  9* 


(41) 


where 

r*m(£)  (^o  +£2  +  2r0£cos^*)1/2.  (42) 


This  control  law,  which  is  referred  to  as  the  constant  in-plane  angle  control  method ,  results 
in  the  existence  of  an  equilibrium  point  of  (33)-(36).  Moreover,  the  associated  equilibrium 
point  of  system  (33)-(36)  will  then  be  (<?*,  0,0,0),  where  9*  is  the  desired  in-plane  angle. 
V.  Stability  Analysis  of  the  TSS  During  Retrieval 
Suppose  for  simplicity  that  i  <  0  throughout  retrieval.  From  Eq.  (41)  we  have 


<  0 


fir0 


(1 


-)  sin  9*  <  0. 


Denote  by  £{  the  initial  (pre-retrieval)  tether  length.  Then  the  condition  for  i  <  0  is  that 
9*  satisfies  either  0  <  9*  <  f  or  — 7r  <  9*  <  9\  (see  Figure  2),  where  9\  —  9\{£j)  is  such 
that 

cos  9\  —  ,  7r  <  9\  <  ~.  (43) 

2r0  2 
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Figure  2.  Retrieval  Regions  for  9*  with  <j>*  =  0 

From  the  discussion  above,  we  have  £  <  0  and  £  >  0  during  retrieval  by  constraints. 
In  addition,  £  =  0  occurs  only  at  £  —  0.  Hence,  £  will  approach  0  asymptotically. 

Denoting  9  :=  9  —  9* ,  the  linearized  system  state  equations  at  the  equilibrium  point 


(9* ,  0, 0, 0)  are  found  to  be 

9  =  we  (44) 

p  ao  2„5  _  p 

ue  =  (2fi-  cot  9 *  +  (f^(£y)5  ^  ~  2~f6  (45) 

<j>  =  w</>  (46) 

Cjj,  =  (— fl2  +  cot  9*)(j>  —  2^0  (47) 


The  linearized  model  of  the  system  (33)-(36)  is  therefore  seen  to  decouple  into  the  subsys¬ 
tems  (44)-(45)  and  (46)-(47).  It  follows  from  Lemma  2  and  Theorem  3  below  that  these 
two  decoupled  subsystems  are  unstable  for  constant  angle  retrieval  for  any  potential  value 
of  6*  (see  Figure  2),  and  that  the  original  system  (33)-(36)  is  then  also  unstable. 
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Lemma  2.  If  there  is  a  constant  8  >  0  such  that  b(t)  >  8  for  all  t  >  to,  then 
unstable  for  the  system 

Ai  =  x2 

x2  =  a(t)x1  +  b(t)x2. 

Proof:  By  the  Abel-Jacobi-Liouville  theorem  (e.g.,  [10]),  we  have 

det(4>(t,  f0))  =  exp(  f  b(s)ds), 

Jt0 

where  $(t,fo)  is  the  state  transition  matrix  of  the  linear  system  (48)-(49).  From  the 
hypothesis,  it  now  follows  that  det  4>(f,to)  approaches  oo  as  t  — »  oo.  Hence,  the  origin  is 
unstable. 

□ 

Theorem  3.  If  r0  >  £i,  then  the  origin  is  unstable  for  the  linearized  system  (44)-(47) 
during  constant  in-plane  angle  retrieval,  and  (6*  ,0,0,0)  is  an  unstable  equilibrium  point  of 
the  original  nonlinear  system  (33)-(36). 

Proof:  First,  we  show  that  there  is  a  constant  e  >  0  such  that  ^  >  e  for  all  t  >  t0.  By 
the  discussion  above,  the  tether  length  i  approaches  0  asymptotically.  In  addition,  t  —  0 
when  £  =  0.  Applying  L’Hopital’s  Rule  to  (41)-(42),  we  obtain 

p  dl  3 

lim  -  =  lim  —  =  —  -Pi  cos  9*  sin  9*  ^  0, 
t—*o  £  £->0  d£  2  r  ’ 

for  any  admissible  9*  (see  Figure  2).  Moreover,  it  is  obvious  that  —  |  >  0,  for  all  t  >  to 

during  retrieval.  Thus,  -j-  is  bounded  and  e  exists. 

Next,  we  check  the  signs  of  the  coefficients  multiplying  ooe  and  in  Eqs.  (45)  and 
(47),  respectively.  Since  0  <  e  <  ^  for  all  t  >  to,  Lemma  2  implies  the  origin  is  unstable 
for  the  two  subsystems  (44)-(45)  and  (46)-(47).  This  implies  instability  of  the  equilibrium 
point  (9*,  0,0,0)  of  system  (33)-(36)  during  constant  in-plane  angle  retrieval. 

a 


the  origin  is 

(48) 

(49) 
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VI.  Stability  Analysis  of  the  TSS  During  Deployment 

In  this  section,  we  consider  application  of  the  constant  in-plane  angle  strategy  of 
Section  IV  to  subsatellite  deployment.  For  simplicity,  suppose  that  t  >  0  for  all  t  >  to. 
Since  £  is  always  positive  in  this  consideration,  one  might  expect  that  the  tether  length 
t  increases  without  bound.  In  reality,  only  a  finite  final  tether  length  is  meaningful  for 
deployment.  Stability  of  the  TSS  is  hence  only  considered  in  a  finite  time  interval,  where 
Liapunov  stability  criteria  cannot  be  employed.  Results  from  finite-time  stability  shall  be 
applied  to  study  the  behavior  of  this  system. 

Basic  definitions  and  conditions  for  finite-time  stability  are  given  in  Section  VI.  1.  Then 
these  finite-time  stability  criteria,  especially  the  contractive  stability  criteria,  are  used  to 
study  the  stability  of  the  TSS  during  constant  in-plane  angle  deployment.  In  addition 
to  the  proof  of  stability  of  deployment,  a  switching  type  control  law  combining  constant 
angle  deployment  and  station-keeping  control  is  proposed  to  achieve  asymptotic  stability. 
Details  of  this  are  given  in  Section  VI. 2. 

VI. 1.  Results  on  Finite-Time  Stability 

Consider  a  system  given  by 

x  =  f(t,x),  (50) 

where  /  :  T  x  Rn  — »  Rn  and  T  :=  [to, to  +  T)  for  some  t0  6  R,  T  >  0.  Let  x0  denote  the 
initial  condition  of  (50)  at  to,  and  let  </>(t;io,#o)  be  the  solution  of  (50)  at  time  t  satisfying 
the  initial  condition.  Then  we  have  the  following  definitions  [14]. 

Definition  1.  System  (50)  is  finite-time  stable  with  respect  to  ||  •  ||),  cn  <  ft,  if  for 

every  trajectory  <^(t;to,xo)  with  ||xo|]  <  a,  we  have  ||<^(t;to,^o)||  <  ft,  V  t  6  F. 

Definition  2.  System  (50)  is  uniformly  finite-time  stable  with  respect  to  T,  |j  •  ||), 

a  <  ft,  if  for  any  trajectory  <f>(t-,s,x)  with  ||x||  <  a,  V  s  £  T,  we  have  \\<p(t]s,x) ||  <  /?,  V 

t  €  r. 

Definition  3.  System  (50)  is  quasi- contr actively  stable  with  respect  to  (a,7,r,  ||  •  |l),  7  <  «, 
if  for  any  trajectory  <j>(t;  t0,  ®o)  with  |[x0||  <  a,  there  exists  ti  G  T  so  that  \\<j>(t;  t0,  a;0)||  <7, 

V  t  €  [ii,  to  +  T). 
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Definition  4-  System  (50)  is  contractively  stable  with  respect  to  (a,  &7,r,  IHI)>  7  <  «  <  P, 
if  it  is  finite-time  stable  with  respect  to  («,/?,  I\  ||  •  ||)  and  quasi- contractively  stable  with 

respect  to  (0,7,  F,  | f  -  ||). 

For  given  a,/3,T,  and  ||  •  ||,  a  necessary  and  sufficient  condition  for  uniform  finite-time 
stability  can  be  stated  as  follows. 

Lemma  3  [14].  System  (50)  is  uniformly  finite-time  stable  with  respect  to  (a,0,T,  ||  •  ||), 
a  <  0,  if  and  only  if  there  exists  a  continuous  function  V(t,x)  such  that 


v(t,x)<  0,  VxeB(p),  ter,  (5i) 

v&(*0  <  v£(h),  v  t2  >  tu  V  s  <  a,  fi,*2  e  r,  (52) 

where 

B{0)  :=  {a;  :  ||rc||  <  /?},  (53) 

||  •  ||  denotes  any  norm  on  Rn,  B{0)  denotes  the  closure  of  B(/3),  and 

Vfiiit)  :=  sup  V(t,x),  (54) 

IMI=« 

<65) 

11x11=0- 


Here,  V(t,x)  is  the  time  derivative  of  V(t,x)  along  trajectories  of  system  (50). 

Stability  properties  of  a  system  may  be  investigated  without  reference  to  the  specific 
bounds  on  the  states  (i.e.  a,  ft  and  7).  In  the  following  lemma  and  theorem,  two  sufficient 
conditions  are  introduced  for  this  type  of  stability.  These  provide  a  means  for  finding  the 
associated  bounds  a,0,7.  Lemma  4  gives  a  sufficient  condition  for  uniform  finite-time 
stability.  Theorem  4  then  gives  a  relationship  among  T,  a,  0  and  7  providing  a  sufficient 
condition  for  contractive  stability. 

Lemma  4.  System  (50)  is  uniformly  finite-time  stable  with  respect  to  (a,/3,T,  ||  •  ||)  for 
any  given  a  and  0  with 


if  there  exist  r  >  0  and  a  continuously  differentiable  function  V(t,x )  with 

V (t,  x)  <  0, 

&ilM|2  <  V(t,x)  <  fc2||a:||2,  ”  (57) 

for  all  x  G  -B(r),  t  G  F.  Here,  0  <  ki  <  and  the  norm  used  is  the  Euclidean  norm. 
Proof.  The  result  follows  directly  from  condition  (52)  of  Lemma  3. 


o 

In  the  next  theorem,  we  introduce  a  condition  on  (50)  and  a  relationship  among  T,a,P 
and  7  guaranteeing  finite-time  contractive  stability. 

Theorem  4.  System  (50)  is  contractively  stable  with  respect  to  (<*,/?, 7, T,  ||  •  ||)  for  any 
triple  a,  /?,7  with 

•exp(-lT)  si,<a<  &  <68) 


if  there  exist  r  >  0  and  a  continuously  differentiable  function  V(t,x )  satisfying  the  condi¬ 
tions 


h\\x\\2  <V(t,x)  <  k2\\x\\2,  (59) 

h\\x\\ 2<-V(t,x),  (60) 


for  all  x  G  -B(r),  t  G  T.  Here,  >  0,  i  —  1,  2,  3,  ||x||  is  the  Euclidean  norm,  and  the  time 
interval  length  T  is  such  that 


T>f-  lnfi 

ki 


(61) 


Proof:  Condition  (60)  implies  that 

V(t,x)<  0,  V  x  G  B(p),  t  G  r. 

Hence,  it  is  implied  by  Lemma  4  that  (50)  is  uniformly  finite-time  stable  with  respect 
to  (a,/?, T,  ||  •  ||)  for  any  a,  P  satisfying  condition  (58).  Next,  we  prove  quasi- contractive 
stability  of  the  system.  From  conditions  (59)  and  (60),  we  have 

V(t,  x)  <  —  ~V(t,x),  VxG  B(r ),  t  G  T. 
ki 
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Hence, 


£  _ 

V(t,cj)(t-,to,xo))  <  V(t(hx0)eyLv(--f-(t -t0)),  V  x0  G  R(r),  t  e  F. 

K2 

Then  it  follows  from  (59)  that 

||^(f;f0,zo)||2  <  ^IM|2exp(-^(t  -  t0)),  V  x0  G  B(r),  t  E  T. 

rCi  K  2 


Thus,  there  exists  a  G  F  so  that  ||^(t;t0,*o)ll  <  7,  V  t  G  [<i ,  #o  +T)  when  conditions  (58) 
and  (61)  hold.  Then  according  to  Definition  3,  system  (50)  is  quasi-contractively  stable 
with  respect  to  (o;,7,r,  ||  •  |j)  for  any  a:, 7  satisfying  (58). 


□ 


VI. 2.  Application  to  Deployment 


In  the  following  discussion,  we  consider  the  deployment  of  the  subsatellite  of  the 
tethered  satellite  system.  For  simplicity,  let  i  >  0  throughout  deployment.  By  Eq.  (41), 
we  have 


t  >  0  4=^ 


f2r0 


(1- 


„3 

r0 


-)sin0*  >  0 


From  the  discussion  above  and  Eq.  (42),  the  condition  on  Q*  for  i  >  0  is  that  either 
62 (if)  <  0*  <  7r,  or  —  j  <  8*  <  0  (see  Figure  3),  where  ^2(^7)  solves 


cos  0 2 


—  ,  0  <  02  <  7T 

2  r0 


(62) 


and  if  is  the  desired  post- deployment  tether  length. 
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Figure  3.  Deployment  Regions  for  6*  with  <j>*  =  0 

Two  strategies  for  deployment  are  considered  here.  The  first  consists  of  the  constant 
in-plane  angle  control  law  for  deployment,  and  the  second  involves  following  the  constant 
in-plane  angle  control  law  followed  by  a  stabilizing  station- keeping  control  once  the  desired 
in-plane  angle  is  close  enough  to  0  radians  or  7r  radians. 

Strategy  1.  Constant  Angle  Control  Only 

We  now  consider  application  of  the  constant  in-plane  angle  control  law  discussed  above 
to  subsatellite  deployment.  In  the  following,  £f  denotes  the  desired  final  tether  length  and 
£;  denotes  the  initial  tether  length  supported  by  a  boom. 

From  (41),  we  have 

£  flrosin#*  n  r0  ^ 

ft  sinfl*  ro  (r^(£))6-r60 

2  £  (r^(^))3  [rg  +(r^(£)f] 

>  —0,  sing*  cos  6*  >  0  (63) 

for  any  9*  €  Sd  and  £j  <£<  ^r0,  where 

Sd  :=  {0*j  —  0.68  radians  <  6*  <  0,  or  2.5  radians  <  6*  <  z  radians}. 
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Hence,  f  is  bounded  below  for  all  U  <  t  <  if  and  #*  £  Sd,  and  similarly  for  L  Thus,  for 
any  if  >  t  will  increase  past  If  at  some  T  >  0.  Theorem  5  below  asserts  that  the  system 
will  be  finite-time  contractively  stable  during  deployment  over  the  interval  [i0,  to  +  T),  near 
the  equilibrium  point  (#*,0,0,0)  with  #*  £  Sd- 

Theorem  5.  Suppose  O  <  1,  lf  <  ^r0,  and  V  :=  [t0,t0  +  T).  There  is  an  r  >  0  such 
that  system  (33)-(36)  is  finite-time  contractively  stable  with  respect  to  0,/?,7,r,||  •  ||)  at 
the  equilibrium  point  (#*,0,0,0)  for  any  a,  /?,  7  and  T  satisfying  (58)  and  (61),  if  either 
of  following  two  conditions  on  the  desired  in-plane  angle  #*  holds: 

(i)  —0.68  radians  <  #*  <  0 
(ii)  2.5  radians  <  #*  <  7r  radians 

Proof:  Let  m  f2sin2#*,  then  it  is  clear  from  (63)  that 


-~j  <m<  0,  V  t  e  r, 

if  either  of  conditions  (i)  and  (ii)  holds  and  ~r0  <  If  <  £.  Invoking  the  finite-time  stability 
criteria  given  in  Section  VI. 1,  the  stability  of  the  TSS  during  constant  angle  deployment 
can  be  proved  as  follows. 

Using  a  general  construction  [18]  for  a  class  of  systems  of  the  form  (48)-(49),  we 
prove  the  finte-time  contractive  stability  of  (33)- (36)  during  deployment  by  employing  the 
Liapunov-like  function 

V(t,#,a>e,0,^)  =  (^  +  — -^-^)#2  +  2 9lo6  -  —cog 

z  m  m 

,2£  ri2 (t).  ,9  ,  1  9 

+(7  +  “^r^  +2M (64) 


where 


P  qo2r5 


sin2  #* 


ri2 (t)  :=  — fl2  +20-  cot#*. 


(65) 

(66) 
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Then  corresponding  to  the  original  system  (33)-(36),  we  have 


2  i 


V(t,  +  n±(t)(f)  +  2(1  -f-  — — )  •  (ui^  +  tu^) 


Vo  , 


+  2(e-=i)/1(0  +  2(^-^)/2(0, 

m  m 


(67) 


where 


h(t)  =  -^ft  +  2tan <f>(u0  +  Q)u4,  -  ~  ~  TTTT^ )  ~  ni(t)d,  (68) 


21. 

£ 


£cos  (j> 


„3 


1  Q?rn 

Mt)  =  —  2  sin(2</»)(u7e  +  ^)2 - ^-cos(0*  +  6f)sin^(l  -  -  n2  (*)</>,  (69) 


and 


.  .  .  .  d  .  2th  1  dn\(t ) 

/  \  „  t  \  d  2£.  1  dn2(t) 

«.(!)- 2n,(f)  +  3i{7}  +  -.-3ii. 


(70) 

(71) 


First,  consider  the  case  in  which  6*  satisfies  condition  (i).  After  some  calculations 
using  (64)-(66),  we  find  that  there  exist  k\t\,  k1>2  >  0  (given  in  the  Appendix)  such  that 

h, llMI2  <  V(i,0,u>tf,^,w*)  <  kl)2\\x\\2 ,  V  t  e  T,  (72) 


where  x  =  (8,ujq,  (f>,to(j,)T  and  the  norm  indicated  is  the  Euclidean  norm.  Moreover,  by 
choosing  k\ ^  0.132122  and 


r  =  su 

«er 


fllll  IX  I  /  mfel,3|k||  _J  IX  I  /  m^l,3|kll1 

p(lwl  :  1/11 5  K^rt)  md  ^2(^1 J>' 


(73) 


we  have 


-I7(*,0,w0,^,u>0)  >  e  r,  xe  B(r). 


(74) 


Thus,  conditions  (59)-(60)  are  satisfied  and  the  conclusion  follows  from  Theorem  4. 
Similarly,  for  the  case  in  which  9*  satisfies  condition  (ii),  we  have 

k2M\2<v(tJ,ve,^v<t>)<k2,2\\x\\\  Vie r,  (75) 
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where  fc2)i,  ^'2,2  >  0  are  also  specified  in  the  Appendix.  By  choosing  &2)3  :=  0.0442fi2  and 

(76) 


. . .  ,  ,,  ,  „  mi2ia||i|| 

r  =  sup{  11*11 :  i/n  <  20ST=rry  and  1/21  s  ^n}’ 


we  guarantee  that 


-V(t,9,uit,<l>,ui<i,)  >  i2,3!l*li2i  Vi  e  r,  xe  J3(r). 


(77) 


The  conclusion  again  follows  from  Theorem  4. 

□ 

The  finite-time  contractive  stable  regions  of  the  desired  in-plane  angle  0*  for  constant 
in-plane  angle  deployment  are  given  in  Theorem  5.  In  addition,  a  relationship  between  the 
time-interval  T,  the  bound  of  initial  disturbances  and  the  final  contracted  region  is  set  up 
in  Theorem  4.  Furthermore,  the  simulation  results  given  in  Section  VII. 2  shows  that  the 
criteria  given  in  Theorem  5  are  not  vacuous. 

Strategy  2.  Station- Keeping  Control  Included 

A  tension  control  law  has  been  designed  in  Section  III  to  regulate  the  tether  length 
at  the  point  where  the  out-of-plane  angle  (f>  =  0  and  the  in-plane  angle  9  —  0  or  {9  —  7r). 
Provided  by  Theorem  5  and  the  control  strategies  given  in  Theorems  1  and  2  (which  design 
tension  control  law  with  linear  state  feedbacks  and  nonlinear  state  feedbacks  from  the  Hopf 
bifurcation  theorem),  a  switching  control  law  for  deployment  can  be  set  up  as  follows: 

Step  1.  Apply  the  constant  angle  control  law  (41)  for  the  first  step  subsatellite 
deployment,  in  which  the  desired  in-plane  angle  9*  satisfies  the 
conditions  of  Theorem  5  and  closes  to  0  radians  (or  7 r  radians). 

Step  2.  Apply  the  tension  control  law  given  in  Theorem  1  (or  Theorem  2) 
once  the  tether  length  is  sufficiently  near  the  desired  length  If. 

It  is  supported  by  Theorem  5  that  the  initial  disturbance  states  of  the  TSS  can  be 
contracted.  Especially,  with  the  desired  in-plane  angle  is  sufficiently  near  0  or  7r  radians, 
the  system  states  of  the  TSS  are  expected  to  be  contracted  within  the  domain  of  attraction 
of  the  station- keeping  control  mode,  after  first  step  constant  angle  deployment.  Hence,  the 
tether  length  should  be  regulated  to  the  desired  length  after  switching  to  the  stabilization 
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control  proposed  in  Section  III  and  [13]  when  the  difference  of  the  tether  length  with  respect 
to  the  desired  value  is  small  enough.  Thus,  the  TSS  is  concluded  to  be  asymptotically  stable 
for  deployment  by  using  this  switching  control  law.  Simulation  results  of  a  typical  system 
given  in  Section  VII.2  illustrate  this  conclusion. 

VII.  Simulation  Results 

Many  simulation  examples  for  tethered  satellite  systems  in  the  station-keeping  mode 
have  been  presented  in  [13].  In  this  section,  we  present  simulation  results  only  for  deploy¬ 
ment  and  retrieval. 

A  typical  TSS  with  following  characteristics  is  considered  : 

•  Orbital  radius  ro  =  6598  km, 

•  Subsatellite  mass  m  =  170  kg, 

•  Orbital  angular  velocity  =  0.0011781  radian/second. 

In  the  following  discussions,  9  =  9  —  6*  denotes  the  differential  of  the  in-plane  angle, 
i f  denotes  the  desired  final  tether  length,  t  =  l  —  if  denotes  the  differential  of  the  tether 
length  and  Ft  denotes  the  applied  tension  control  force. 

VII.  1.  Retrieval 

As  discussed  in  Section  V,  the  set  of  candidate  in-plane  angles  for  constant  angle 
retrieval  is  given  by 

Sr  :=  {0|O  <  9  <  —  or  —  n  <  9  <  (£;)}, 

where  #i(£j)  is  defined  in  (43).  Let  the  initial  disturbances  of  the  system  be  (f>  =  0.01 
radians,  9  =  —0.01  radians,  and  u>q  =  =  0.  The  initial  tether  length  is  assumed  to 

be  10  km.  It  is  observed  from  Figures  4  and  5  that  the  equilibrium  point  (9*,  0,0,0)  is 
unstable  during  retrieval  with  the  desired  in-plane  angle  9*  =  —3.0  radians  and  9*  =  —1.6 
radians,  respectively.  As  mentioned  in  [13],  in  reality,  since  the  tether  is  not  I'eally  rigid, 
the  applied  tension  control  force  can  not  be  positive  to  rule  out  compression.  However, 
Figure  5  (d)  shows  that  a  positive  tension  control  force  Ft  occurs  in  some  time  interval. 
Thus,  if  a  constant  angle  control  law  is  applied  during  retrieval,  then  not  only  will  the 
system  be  unstable,  but  also  tether  compression  may  occur. 
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It  is  also  found  that  whenever  the  desired  in-plane  angle  9*  £  Si  :=  {—2.1  radians 
<  9*  <  9i(£i)}  for  constant  angle  retrieval,  the  applied  tension  control  force  Fi  can  have 
positive  value  in  some  time-interval,  i.e., compression  may  occur.  The  system  response  of 
9*  =  —2.1  radians  for  constant  retrieval  is  depicted  in  Figure  6,  where  Fi  is  found  (see 
Figure  6  (d))  to  be  very  close  to  0  at  some  time  instant  but  never  greater  than  0. 

Similar  simulation  results  are  found  for  the  region  0  <  9*  <  j  for  constant  angle 
retrieval.  The  equilibrium  point  (9*,  0,0,0)  is  found  unstable  during  retrieval  and  the 

compression  of  the  tether  may  occur  in  case  1.0  <  9*  <  j.  The  system  responses  are  not 

shown. 

VII. 2.  Deployment 

According  to  Theorem  5,  the  set  of  candidate  9*  for  stable  deployment  is 

Sd,  =  {#|  —  0.68  radians  <  9  <  0,  or  2.5  radians  <  9  <  ix  radians}. 

Let  the  initial  disturbances  of  the  system  be  <f)  =  0.01  radians,  9  =  —0.01  radians,  and 
loq  ~  u)(j>  —  0.  The  initial  tether  length  is  assumed  to  be  l\  —  10  m,  which  is  provided  by 
a  boom.  First,  the  system  response  during  deployment  (applying  constant  in-plane  angle 
control  only)  are  depicted  in  Figures  7  and  8,  with  9*  —  —0.68  radians,  and  9*  =  2.5 
radians,  respectively.  It  is  observed  from  the  system  responses  that,  for  example,  the 
differential  of  the  in-plane  angle  9  and  the  out-of-plane  angle  (j),  decay  during  deployment. 

The  switching  control  strategy,  which  involves  both  constant  angle  control  and  station¬ 
keeping  control,  is  applied  to  deploy  a  subsatellite  from  the  satellite  with  the  desired  final 
tether  length  if  —  10  km.  The  first  example  considers  to  deploy  the  subsatellite  upward 
(i.e.,  away  from  the  Earth)  by  applying  constant  angle  control  with  9*  —  —0.015  for  first 
260,500  seconds,  and  applying  the  station-keeping  control  thereafter.  The  applied  tension 
control  force  for  station-keeping  is  governed  by 

Ft  =  - m(U  +  M  +  h2l),  (78) 

where  U  —  0.041019,  hi  =  3.1fi2  and  h2  =  0.0034.  The  responses  of  the  system  during 
constant  angle  deployment  are  shown  in  Figure  9.  At  time  t  =  260, 500  seconds,  we  have 

•  the  out-of-plane  angle  (f>  =  —7.01636  x  10-6  radians  and 
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<f>  —  1.70633  x  10-8  radian/second 

•  the  in-plane  angle  9  —  —0.0150051  radians  and  9  =  5.61812  x  10~10  radian/second 

•  the  actual  tether  length  i  =  9.97617  km  and  £  =  2.63603  x  10-4  km/second. 

With  these  values,  the  applied  tension  control  law  is  switched  to  the  station-keeping 
control  and  governed  by  Eq.  (78).  The  system  responses  governed  by  (78)  are  depicted  in 
Figure  10. 

Another  example  for  deploying  subsatellite  downward  (i.e.,  toward  the  Earth)  is  im¬ 
plemented  by  applying  constant  angle  control  for  the  first  235,300  seconds  with  6*  =  3.125 
radians,  then  switched  to  the  station-keeping  control  governed  by  Eq.  (78).  At  time 
t  =  235, 300  seconds,  we  have 

•  the  out-of-plane  angle  <f>  =  —2.01378  x  10-6  radians  and 
<j)  =  —2.35517  x  10~8  radian/second 

•  the  in-plane  angle  9  —  3.1250  radians  and  9  =  8.96566  x  10~9  radian/second 

•  the  actual  tether  length  i  =  9.88531  km  and  £  =  2.90670  x  10~4  km/second. 

The  system  responses  are  shown  in  Figures  11  and  12. 
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Appendix 

The  values  of  kij,  i  —  1,2  and  j  =  1,2  are  as  given  below. 

0.295575407  esc2  20* 


.  _  h  +  h  +  yj (h  +  h)2  +  4 

1,2  _  2  ’ 

,  0.498328311  esc2  29* 

"■2,1  =  - 7 - > 

k2,2 

?4  +  ^5  +  y/(h  +  h)2  +  4 

k2,2  =  - - - , 


where 

h 
h 

h 
h 
h  — 

h  = 


—30  cos  6*  sin  6*  —  0.50  sin  9* , 

02(3  cos2  9*  +  0.5  cos  (9*)  +  1 
Osin  29* 

02(3  cos2  9*  +  0.5  cos  9*)  -  1 
O  sin  29* 

—3.689610  cos  9*  sin  9* , 

3.6896061O2  cos2  9*  +  1 
O  sin  29*  ’ 

3-6896061O2  cos2  9*  -  1 
O  sin  29* 
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Figure  5.  Simulation  Results  for  Constant  Angle  Retrieval 
with  9*  =  —1.6  radians 
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Figure  6.  Simulation  Results  for  Constant  Angle  Retrieval 
with  9*  —  —2.1  radians 


Figure  7.  Simulation  Results  for  Constant  Angle  Deployment 
with  9*  =  —0.68  radians 
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Figure  8.  Simulation  Results  for  Constant  Angle  Deployment 
with  6*  =  2.5  radians 
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Figure  9.  Simulation  Results  for  Constant  Angle  Deployment 
with  6*  =  -0.015  radians 
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Figure  12.  Simulation  Results  for  Station-Keeping 
with  9*  =  7r  radians 
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